Multiple solutions of coupled-cluster equations for PPP model of [10]annulene 
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Multiple (real) solutions of the CC equations (corresponding to the CCD, ACP and ACPQ meth- 
ods) are studied for the PPP model of [10]annulene, CioHio. The long-range electrostatic inter- 
actions are represented either by the Mataga-Nishimoto potential, or Pople's R~ l potential. The 
multiple solutions are obtained in a quasi-random manner, by generating a pool of starting ampli- 
tudes and applying a standard CC iterative procedure combined with Pulay's DIIS method. Several 
unexpected features of these solutions are uncovered, including the switching between two CCD 
solutions when moving between the weakly and strongly correlated regime of the PPP model with 
Pople's potential. 



I. INTRODUCTION 

The coupled-cluster (CC) method is one of the ba- 
sic tools for studying the electronic-correlation effects in 
many-electron systems. Due to their nonlinearity, the CC 
equations may have multiple solutions; this problem was 
first analyzed by Zivkovic and Monkhorst [Q. Recently, 
Kowalski and Jankowski [|| applied to the CC equa- 
tions a new powerful mathematical tool the homotopy 
(continuation) method Q. By using this tool they were 
able to find the complete sets of solutions for the CCD 
and the CCSD equations corresponding to some four- 
electron four-orbital systems, known as the H4 and P4 
models. Very recently, Podeszwa et al. || used the homo- 
topy method in a study of multiple solutions to the CCD 
equations for the Pariser-Parr-Pople (PPP) model |(| [t[] 
of benzene (6 electrons, 6 orbitals). The available data 
on the CC multiple-solution problem is still very limited. 
The purpose of this paper is to get some insight into 
the multiple solutions of the CCD equations (as well as 
the equations of some modified CCD methods) for the 
PPP model of [10]annulene (10 electrons, 10 orbitals), 
a system which is already too big to be treated by the 
homotopy method. 

Benzene and [10]annulcnc arc the first two members 
of the family of cyclic polyenes (annulenes) CmHm, 
where M = 4mo + 2, mo = 1,2,..., and the assumed 
point-symmetry group is Dj\/h- The problem of elec- 
tronic correlations in annulenes, which may serve as 
prototypes of one-dimensional metal, have been thor- 
oughly studied within the PPP model by Paldus and 
co-workers p[ p|, |Ifj||. Due to the high symmetry of an- 
nulenes, the contributions from the single excitations are 
absent in the PPP model, and the CCSD method re- 
duces to the CCD one. Paldus et al. discovered and doc- 
umented the breakdown of the CCD method for these 
systems. They also showed that this breakdown may 
be attributed to the neglect in the CCD method the 
terms providing the coupling between the D and Q exci- 
tations. Paldus et al. devised an approximate coupled- 
pair method, called ACPQ ^ M, in which these coupling 
terms (approximately) canceled certain quadratic terms 
in the CCD equations. Unlike the CCD method, the 
ACPQ method, and a related ACP-D45 method (ACP in 



short) introduced earlier pd| , were shown to perform well 
for the annulenes, being convergent and giving the corre- 
lation energies close to the full configuration-interaction 
(FCI) results. However, a recent CC study [Q of annu- 
lenes, taking into account the double (D), triple (T), and 
quadruple (Q) excitations in the CC operator, showed 
that even the CCDTQ method breaks down for these 
systems when the correlation effects become sufficiently 
strong. Moreover, the ACP and ACPQ methods were 
shown to fail for very large annulenes |l3|]. Thus, the 
annulenes remain a challenging problem for the many- 
electron theories. 



II. FINDING MULTIPLE SOLUTIONS OF CCD, 
ACP, AND ACPQ EQUATIONS FOR 
[10] ANNULENE 

In the CCD method, the CC operator corresponding 
to the double excitations depends on some linear param- 
eters, hereafter referred to as the t 2 amplitudes. The 
CCD equations comprise a set of quadratic equations for 
the unknown ti amplitudes. In the nonorthogonally spin- 
adapted formalism (see, e.g., Ref. |l4|), the number of t?, 
amplitudes is equal to K = n n u (n n u + l)/2, where 
n and n u is the number of occupied and unoccupied 
orbitals, respectively (no symmetry other then the spin 
one is assumed) . According to the Bezout theorem, such 
a set of quadratic equations may have up to 2 K solu- 
tions, a truly astronomical number in most applications. 
Therefore, finding the complete set of CCD solutions 
may be attempted only when the number of t% ampli- 
tudes is very small (e.g., due to the spatial symmetry of 
the many-electron system). In the homotopy study of 
the H4 model || (corresponding to 6 spin- and spatial- 
symmetry-adapted ti amplitudes), the complete set of 
CCD solutions numbered 12, which is distinctly smaller 
than the Bezout limit (64). On the other hand, the ap- 
plication of the homotopy method to the PPP model of 
benzene § (having 8 t 2 amplitudes adapted to spin, spa- 
tial, and time-reversal symmetry) has brought in some 
quite disturbing findings: not only the total number of 
solutions (230) approached the theoretical limit (256), 
but some solutions led to the correlation energies much 



2 



more negative than that of the "standard" CCD solu- 
tion (which approximate the ground-state FCI solution) ; 
in addition all the "nonstandard" solutions were hardly 
similar to the FCI ones. 

In the PPP model of [10]annulcnc, the CCD equa- 
tions adapted to spin, spatial, and time-reversal sym- 
metry correspond to 29 t% amplitudes. The ACP and 
ACPQ equations are obtained by neglecting (ACP and 
ACPQ) and modifying (ACPQ) some of the quadratic 
terms in the CCD equations, see, e.g., Ref. |l3|. In the 
present paper we use the PPP model of annulcnes de- 
scribed in Ref. [T^ and follow the notation used there: (3 
(< 0) eV is the resonance integral of the PPP model, 
and 7(i?) is a potential describing the long-range elec- 
trostatic interactions. We shall report the results ob- 
tained for two variants of the PPP model: the PPP-MN 
model, employing the Mataga-Nishimoto potential (lj| 
( 7 (i?) = e 2 [i?+e 2 (7 )- 1 ]- 1 ), with 7° = 10.84 eV, and the 
PPP-P model, using Pople's potential @ (j(R) = e 2 /R 
for R ± 0), with 7 ° = 14.25878 eV (e is the elementary 
charge) . 

If the CCD results for benzene || may be extrapolated, 
the CCD (and ACP and ACPQ) equations for [10]annu- 
lene may have more than 10 8 solutions; that makes the 
application of the homotopy method hopeless in this case. 
On the other hand, it was found in Ref. |5| that some of 
the real CCD solutions for benzene were stable within 
their close neighborhoods: such a solution (subject to 
the numerical rounding errors) led to a convergent iter- 
ation process when the direct-inversion-of-the-iterative- 
subspace (DIIS) method Q was used. Thus, one may 
expect that the DIIS method is able to trace various so- 
lutions, if only the starting points are suitably chosen. 

A standard method of solving the CCD equations is 
based on the iterative procedure represented in Eq. (14) 
of Ref. @. By starting from the i2 amplitudes equal 
to zero, one recovers in the first iteration the t 2 am- 
plitudes corresponding to the first order of the M0ller- 
Plesset (MP) perturbation theory; in the subsequent it- 
erations higher-order MP contributions are added. Such 
an iteration procedure, if convergent, furnishes the stan- 
dard CCD solution corresponding (as it is often a priori 
assumed) to the ground state of the system. The CCD 
electronic-correlation energy is then calculated as a cer- 
tain linear function of the converged £2 amplitudes (the 
tij amplitudes correspond to the second-order MP cor- 
relation energy and will be, after Ref. |l2|, referred to as 
the MP2 amplitudes). The above treatment applies also 
to the ACP and ACPQ equations. In annulenes one may 
also use a kind of analytical-continuation procedure to 
find the standard solution for those values of parame- 
ter (3 for which the standard procedure (starting from 
the MP2 amplitudes) does not converge: one uses as the 
starting point the converged ti amplitudes found for the 
sufficiently close (3 value. 

The multiple solutions of the CCD, ACP, and ACPQ 
equations presented in this paper were found in a quasi- 
random manner, by generating a pool of the initial ti 



amplitudes and then applying the DIIS procedure. The 
starting pool included the results of several hundred con- 
verged CCD, ACP, and ACPQ calculations obtained for 
different (3 values and for different forms of function 7(i?) 
(these sets of ti amplitudes indeed looked random). Most 
of the calculations starting from this quasi-random am- 
plitudes either diverge, or converge to the standard so- 
lution, but in some cases a new solution is obtained. 
Certainly, this is not a very efficient method of finding 
the multiple solutions, but since a single CCD calcula- 
tion for [10]annulene usually takes only several seconds, 
a lot of data may be generated that way. We have found 
that the strongly correlated regime (f3 > —1.5 eV) is 
rich in multiple solutions, while in the weakly correlated 
regime (/3 < —3.0 eV) practically only the standard so- 
lution could be reached that way. Therefore, we per- 
formed calculations for several values of p from interval 
(— 1.0 eV, eV). These results were then extended as 
far as possible by a careful application of our analytical- 
continuation procedure. The obtained results are not 
supposed to provide the complete set of real solutions 
for [10]annulene. Nevertheless, they provide a rather in- 
triguing glimpse at the "tip of the iceberg" of the multiple 
solution problem for this system. Let us note that the 
multiple solutions of the ACP and ACPQ equations have 
never been studied before. 



III. RESULTS 

In Figs. |l|-^ we present the electronic-correlation ener- 
gies corresponding to the multiple solutions of the CCD, 
ACP, and ACPQ equations, respectively. These results 
are obtained within the PPP-MN model for —5.0 eV < 
(3 < eV, each solution is represented by a continuous 
line. The line endpoints (not touching the graph's bound- 
aries) indicate the limiting (3 values corresponding to our 
analytical-continuation procedure. It is seen that some of 
the new solutions exist for a very broad range of (3, while 
some can be continued only within a very narrow inter- 
val. The ground-state FCI correlation energies, shown 
as the broken line, are provided for a reference. The so- 
lution line that approaches the FCI line in the weakly 
correlated regime corresponds to the standard solution 
discussed in the previous papers [|[ [l2] . 

The CCD results in Fig. are rather simple: The 
standard-solution line is the only extending for all the 
(3 values. (When the standard iteration procedure is 
combined with the DIIS technique, the standard CCD 
solution for [10]annulene can be obtained for all nega- 
tive values of /3, as shown previously by Paldus et al. M 
by using the Newton-Raphson method.) In the strongly 
correlated regime ((3 > —1.5 eV), the standard CCD solu- 
tion provides a rather poor approximation to the ground 
state, and the corresponding correlation energy falls well 
below the FCI line. Other solutions have the energies 
higher than the FCI result, some of them may represent 
excited states. Most of them can be found only in the 
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proximity of (3 = 0, where the quasidegeneracy effects 
are the strongest. 

The ACP and ACPQ results (Figs. | and | respec- 
tively) display a different and more complicated behav- 
ior. Unlike the CCD standard solution, the ACP and 
ACPQ standard solutions stay close to the FC1 line also 
in the strongly correlated regime (although the corre- 



sponding t2 amplitudes show increasing deviations from 
the FC1 t 2 amplitudes as ->■ @). The ACP and 
ACPQ standard solutions can be found for (3 < by 
applying the standard iteration procedure starting from 
the MP2 amplitudes 0. Quite a few nonstandard so- 
lutions have been found, some of them extending for 
the whole range of [3. Quite unusual is the presence 
of the solution lines lying below the standard-solution 
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TABLE I: PPP-P model of [10]annulene, CCD equations. Correlation energies for two solutions, CCD1 and CCD2; ground- 
state MP2 and FCI results shown for comparison (all energies in eV) . 9 (in deg) and r\ gauge the similarity of a given solution 
to the FCI solution (see text). 

P = -2.500 P = -1.175 P = -0.500 
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line (the "underground" lines). Such exotic solutions 
have been observed in the study of benzene ||; however, 
they were unstable in the DIIS iterative procedure (could 
be found only by the homotopy method). The "under- 
ground" solutions obviously have no physical meaning, 
and it is a little disturbing that they may appear while 
iteratively solving the ACP or ACPQ equations. One 
of the ACPQ "underground-solution" lines is especially 
peculiar: it starts with the lowest energy for 8 = and 
rises quite steeply with the decreasing 8, crossing sev- 
eral solution lines, including the standard-solution line 
(apparently there is no "non-crossing rule" for different 
solution lines). 

Quite intriguing are the solutions of the CCD equations 
corresponding to the PPP-P model, displayed in Fig. || 
together with the ground-state FCI results. We were 
able to find only two solutions, hereafter referred to as 
CCD1 (the standard solution) and CCD2 (the nonstan- 
dard one). Despite the use of the extensive quasi-random 
search, the CCD iterations invariably led to CCD1 or 
CCD2 (these solutions were so strong attractors that the 
iterations rarely diverged). In contrast to the CCD re- 
sults for the PPP-MN model, the iterations starting di- 
rectly from the MP2 amplitudes converge to the CCD1 
solution in the weakly correlated regime, while in the 
strongly correlated regime the CCD2 solution is obtained. 
In Fig. [I] we mark by the squares the endpoints of the 
respective direct-convergence domains. One may thus 
obtain the CCD1 and CCD2 solutions in a determinis- 
tic way, by starting the calculations within these regions. 
Moreover, by applying our analytical-continuation proce- 
dure, we are able to extend a given solution line outside 
its domain of (direct) convergence. 

The CCDl-solution line resembles closely the line cor- 
responding to the standard CCD solution for the PPP- 
MN model, see Fig. [j], while the CCD2 line has no coun- 
terpart in that figure. The CCD1 and CCD2 lines cross 
for 3 = —1.17479 eV, in the vicinity of the FCI-solution 
line. In Table § we present the CCD1, CCD2, MP2 and 
FCI correlation energies calculated for three representa- 
tive (3 values, including the intersection coordinate. In 
addition, we show there the values of parameters 6 and rj 
which provide a measure of similarity between the vector 
of ti amplitudes corresponding to a given solution and 
the vector of ti amplitudes extracted from the linear co- 



efficients of the ground-state FCI wavefunction: 9 is the 
angle between the vectors, and r\ is the ratio of the vector 
lengths (for details, see Ref. |l2| ). 

From Fig. [| and Table | a consistent picture emerges: 
The CCD1 solution provides an acceptable approxima- 
tion to the ground-state FCI solution up to the vicin- 
ity of the CCD1/CCD2 crossing; then, in the strongly- 
correlated region, it deteriorates and hardly deserves the 
name of the standard solution. The CCD2 solution, on 
the other hand, starts as the "underground" in the weak 
correlation regime (bearing no resemblance to the FCI so- 
lution) , but it improves in the strongly-correlated regime, 
and, finally, becomes a slightly better (but still poor) 
approximation to the FCI solution (the CCD2 correla- 
tion energies are surprisingly similar to the FCI ones in 
the strongly correlated regime). Thus, it looks as CCD1 
and CCD2 swap their roles as the standard and "un- 
derground" solutions. The behavior of the MP2 results 
agrees with the observed convergence pattern for itera- 
tions starting from the MP2 amplitudes. Interestingly, 
if one performed a CCD study of the PPP-P model of 
[10]annulene by applying a standard iteration procedure 
(starting with the MP2 amplitudes) at some representa- 
tive (3 values (say, —5.0, —2.5, —1.5, and —0.5 eV), the 
conspicuous switch between the CCD1 and CCD2 solu- 
tions would have been passed unnoticed. By looking at 
the CCD correlation energies alone, one would then pro- 
claimed a very good performance of the CCD method up 
to very small 8 values, contrary to the results found by 
Paldus et al. || for the PPP-MN model. 

The study of the ACP and ACPQ equations in the 
PPP-P model brings in some additional unexpected find- 
ings: the standard ACP solution cannot be found for (3 
between -0.80 and -0.4 eV, while the standard ACPQ 
solution encounters problems already for 8 > —2.53 eV. 
Apparently, in the strongly correlated regime of the PPP- 
P model both methods become unstable. 

Our search of the multiple real solutions of the CCD, 
ACP, and ACPQ equations has brought several unex- 
pected findings: (i) while solving the CC equations by 
iterations, one may arrive at an "underground" solution, 
(ii) in the PPP model, the character of multiple solutions 
is sensitive to the form of the long-range potential (the 
7 function), (iii) different solutions may play the role of 
the standard solution when the strength of the electronic- 
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correlation effects is varied (this may apply to the changes 
in molecular geometry, e.g., to the bond breaking). In 
conclusion: it seems that the problem of multiple solu- 
tions of the CC equations deserves still more attention. 
Some variants of our random-search approach may be 



tried also on larger many-electron systems. 
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